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The dynamics of quasi-particles in repulsive Bose condensates in a harmonic trap is studied in 
the classical limit. In isotropic traps the classical motion is integrable and separable in spherical 
coordinates. In anisotropic traps the classical dynamics is found, in general, to be nonintegrable. For 
quasi-particle energies E much smaller than the chemical potential n, besides the conserved quasi- 
particle energy, we identify two additional nearly conserved phase-space functions. These render 
the dynamics inside the condensate (collective dynamics) integrable asymptotically for E/fi —> 0. 
However, there coexists at the same energy a dynamics confined to the surface of the condensate, 
which is governed by a classical Hartree-Fock Hamiltonian. We find that also this dynamics becomes 
integrable for E/ jj, because of the appearance of an adiabatic invariant. For E/fi of order 1 
a large portion of the phase-space supports chaotic motion, both, for the Bogoliubov Hamiltonian 
and its Hartree-Fock approximant. To exemplify this we exhibit Poincare surface of sections for 
harmonic traps with the cylindrical symmetry and anisotropy found in TOP traps. For E/fj, 3> 1 
the dynamics is again governed by the Hartree-Fock Hamiltonian. In the case with cylindrical 
symmetry it becomes quasi-integrable because the remaining small chaotic components in phase 
space are tightly confined by tori. 
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I. INTRODUCTION 

The experimental realization of Bose-condensates of atoms harmonicaUy bound in magnetic traps call for a 
space-dependent version of Bogohubov's theory, or some modification thereof. Such a theory proceeds by splitting 
the field operator ip{x) and its adjoint in a C-number part ipoi^) Si,nd a residual operator ip{x), 

ip{x) = il}o{x)-\-'f{x) (1) 

and an accompanying decomposition of the Hamiltonian in terms of 0, 1, 2, 3, 4 order in ip, ip^ . The term of 1 order 
in ip, ip~^ is made to vanish by choosing tpQ^x) to satisfy the time-independent Gross-Pitaevskii equation which at 
low temperatures, takes the form 

- 7^VVo(aj) + {U{x) - fi)Mx) + Vo\Mx)\^Mx) = , (2) 
2m 

with the normalization J \ipo\'^d'^x — Nq . Here 

U{x)^'^(ulx'+uy+u.lz') (3) 
is the generally anisotropic harmonic trap potential, 
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Vo = (4) 

TO 

is the strength of the pseudo-potential replacing the true two-particle potential at low energies, with the s-wave 
scattering length a, which is here assumed to be positive. 

For (NQa/do) ^ 1, where do = \/h/muj, uj = (o^jjWyWz)^/^, the solution to the Gross-Pitaevskii equation can be 
determined in the Thomas-Fermi approximation [p| by neglecting the kinetic-energy term 

^^^^^^i^^UixlQ^^^^^^^^^^^ (5) 

In the following we shall choose ■00 as real and positive. The chemical potential is determined from the normaliza- 
tion. The next step is the diagonalization of that part of iJ, which is a quadratic form in ip, (y5+, by a Bogoliubov 
transformation to quasi-particles 

(^(a;) = ^([/,(a;)<i,-l-;(a^)a+) (6) 
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with 

j <ex{\U,{x)\'~\V,{xr)^l (7) 

and 

[a,,aj,] = [a1j.a+] 

[aj,a^,]= 5jj, . (8) 
The second-order part of H is diagonalized by this transformation, if Uj{x) and Vj{x) satisfy the Bogoliubov equations 

Fhp -K{x) \ (U,\_^ ( U, 



with the Hartree-Fock Hamiltonian 



t2 

H^^^-—^^ + U{x)-ii + 2Vo\M^)\'' (10) 
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and the coupling term 



K{x) = Vo\Mx)\'' (11) 

between the two components Uj{x), Vj{x) of a quasi-particle wave-function. Because of the different signs of the Uj, 
Vj components on the right-hand side, they play the role of particle and anti-particle components of the complete 
wave- function. As the equations are symmetric under the particle-antiparticle transformation —E^, Uj — > V* , 

Vj — !■ U* we may define Ej to be non-negative without restriction of generality. Various numerical ||7| and approximate 
analytical js) treatments of these equations are available in the literature. 

In the present paper we wish to study the classical limit of the center-of-mass motion of the quasi-particles. In order 
to discuss the dynamics rather than the eigenstates of the quasi-particles, it is useful to introduce time-dependent 
wave functions via 

with arbitrary coefficients Cj . They satisfy the time-dependent Schrodinger equation 
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For large energies Ej, Ej ^ /z, the classical motion can be interpreted as the center of mass-motion of quasi- 
particle wave packets. For small energies Ej, Ej <^ fj,, such a straightforward physical interpretation of the classical 
quasi-particle dynamics is no longer possible. However, even in this regime, there is still a close mathematical 
relation between the classical and the quantum dynamics, as the classical trajectories are the characteristics of the 
quantum mechanical wave-equation. This is made explicit by the derivation of the classical dynamics as a limit of the 
Schrodinger equation via the Hamilton-Jacobi equation. The Hamilton-Jacobi equation corresponding to eq. ( |l3| ) is 
obtained by the asymptotic ansatz for ?i — > 



U{x,t) \ / ao{x,t) + Q{n) \ ,s{x.t)/h 

v{x,t) ) - {bo{x,t) + o{h) 

with J (Px{\ao{x,t)\'^ — |6o(a^,t)P) = 1. It reduces eq. (|l3|) to the form, to zeroth order in h, 

eHF + 7§- —K \ f ao 



(14) 



Here 



e^^ = ^ + U{x)~ti + 2Vo\Mx)\^ ■ (16) 
Zm 

We may restrict to —E = ^ < in accordance with our restriction on E. To first order in h we obtain 

I ( X ) - 1^- ((^^' ( : )) + i^^-^ ( : ) 4 ( f )(?;)■ 

Here (^^) are the O(?i)-components of the amplitudes in (p^. These will exist, and the expansion will be well-defined, 
only if the left-hand side of is orthogonal on the kernel (^") of the matrix in (^5|), which also appears on the 
right-hand side of ([l7|). This condition gives rise to the conservation law 

^(|aop - |6on + ^V.((|aop + \b„\')VS) = (18) 
which ensures that the normalization condition 

d^xdaoP - IfeoH = 1 



is consistent with the classical dynamics and represents the classical limit of the continuity equation following from 
(O) The zeroth order equation has a nontrivial solution only if the determinant condition 



(19) 

is satisfied, which, observing our sign convention for dS/dt, gives the time-dependent Hamilton-Jacobi equation 

9Six,t) ,^f^dS\^ (20) 



dt V dx 

with the classical Hamiltonian 



Hix, p) = ^eUx,p)-K{xr . (21) 
The time-independent Hamilton Jacobi equation results from the separation 

S{x,t) ^ S{x) - Et (22) 

and reads 

H{x,^)=E. (23) 
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If we can neglect the 0(7i)-corrections ai,6i in ( [l7| ) we obtain from the first order equation separate conservation laws 
for the quasi-particle and anti-quasi-particle densities 

^^\ao{x,t)f + ^V- {\ao{x,t)fVS) =0 

^^\bo{x,t)f~^V.{Mx,t)fVS) =0. (24) 

The classical anti-particle and particle dynamics are therefore just the time-reversed of each other, and the densities 
of both components are separately conserved. 

In the following sections we analyze the classical dynamics described by the Hamiltonian (pT|). 



II. CLASSICAL QUASI-PARTICLE DYNAMICS 



For the case of isotropic harmonic traps angular momentum is conserved and the quasi-particle dynamics is inte- 
grable and separable in spherical coordinates. This case is discussed in where it is made the basis of a semiclassical 
quantization procedure. Therefore, in the following we concentrate on the analysis of the case of anisotropic harmonic 
traps in the limit where the Thomas-Fermi approximation applies. In the present section we shall assume cylindrical 
symmetry of the trap 

U{x)^^ix' + y^) + ^z' (25) 

In the experiment |l| > ojo, namely, (w^/wq)^ ~ 8. As the parameter denoting the anisotropy of the potential 
we introduce e by e — 1 — (wq/^z)^ i which is the numerical excentricity of the Thomas-Fermi surface n = U{x), a 
rotational symmetric ellipsoid. This two-dimensional surface is the boundary of the condensate. 

Our problem has a characteristic energy, namely the chemical potential. Thus, the second relevant parameter of 
the classical motion is the ratio E/ii. We note that measuring the energy in units of /z, coordinates, momenta and 
time in units of 



'^o = \ 2 ' Po^ V "^mfi , to^ujQ^ (26) 

respectively, the dimensionless Hamiltonian can be put in a form, which depends only on the anisotropy parameter e. 
This shows that condensates with the same anisotropy but with different chemical potential behave similarly in the 
classical description, if the physical quantities are scaled appropriately. 

In the isotropic case 6 = the classical dynamics are completely integrable. As three independent constants of 
motion we can choose the energy E , the modulus of the angular momentum and its z-component Lz ■ As we keep 
rotational symmetry around the z-axis in the anisotropic case e 7^ the Lz and of course the energy are still conserved 
quantities, whereas the total angular momentum considered here is no longer a constant of motion. Thus, in the 
following we shall investigate the classical behaviour of this three degrees of freedom system depending on the two 
constants of motion E and Lz, and we address the question wether the dynamics are integrable or chaotic. 

Let us introduce the usual cylindrical coordinates p — \J +2/^, z and Because of the rotational symmetry 
around of the z-axes the angle is a cyclic variable. In cylindrical variables the Hamiltonian has merely two degrees 
of freedom p and z, Lz just enters as a parameter. Certain conditions have to be satisfied as can be seen from the 
Hamiltonian in the region outside the condensate. For E > p the condition E + p > ujqLz has to be guaranteed, for 
E < p we must have E > (ujoLz)'^ /Ap . 

The dynamics of this two dimensional system we can visualize by Poincare cuts, see Figj^. For different energies 
we observe different dynamical behaviour. For E > (oj^LzY /'^p > p two different kinds of trajectories can occur 
typically. If the repulsive effective potential in p-direction due to the angular momentum Lz is strong enough, the 
particle cannot enter the condensate and is only moving in the harmonic potential of the external trap. The motion in 
an anisotropic harmonic potential is completely integrable, as a third constant of motion we can choose the energy in 
the z-degree of freedom Ez = pl/2'm + mLOzZ^ . These trajectories, which are not perturbed by the condensate, can 
be seen as the integrable tori around the fixed point of the Poincare map in the centre of Fig.^, which is the periodic 
orbit moving only in z and (p directions. If the particle enters the condensate, Ez is no more a conserved quantity. 
Nevertheless for energies large compared to the chemical potential also those trajectories are still quite similar to 
unperturbed motion. Typically the trajectories are confined to thin stochastic layers separated by each other by 
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integrable tori. No Arnold diffusion occurs, as usually for a system of two, not of three degrees of freedom. At high 
energies the system behaves quasi integrable. The influence of the condensate can be taken as a small perturbation 
to the motion in the external potential. 

For energies in the range 10 > E/ ^ > 0.1 (Fig.|l|b) we typically observe a mixed phase space. The fixed point is 
now inside the condensate, but does not loose its stability. The detailed structure depends on the parameters chosen. 
Already for small anisotropy (e^ = 0.2) a relevant part of phase space can be chaotic. This shows that for energies 
of the order of the chemical potential the isotropic case with its integrable dynamics is an exceptional rather than a 
typical situation. If < ^ all trajectories move inside and outside the condensate. 

For energies small compared to the chemical potential E < O.l/u (Fig.|l|c) the chaotic part of phase space decreases 
again and is restricted to a thin layer separating and surrounding two regular islands, corresponding to two stable 
fixed points separated by an unstable one. Most orbits seem to lie on integrable tori. This suggests that the system 
has an integrable regime in the limit of small energies. 

This limit corresponds to the hydrodynamical regime [|lO| investigated in several contexts. In a bulk case, when 
there is no external potential U (x) the lowest lying excitations are phonons with linear wave-number dependence. 

Numerically we have found that tending with the energy to zero, keeping /i fixed the range of the classical motion 
outside the condensate for trajectories starting inside is getting smaller and smaller and in the limit the motion 
is confined to the region inside of the Thomas-Fermi surface. Starting trajectories from the same point inside the 
condensate under the same direction and changing only the modulus of Cartesian-momentum we have found that 
they differ from each other only in a thin region near the boundary whose width scales with the energy. Lowering the 
modulus of the initial momentum to zero they tend to a well-defined limiting trajectory. This can be clearly seen in 
Fig.^. In the isotropic case this is the limit E/ h^loqL/ ^ 0, keeping the ratio L/E fixed. In the following section 
this 'hydrodynamic regime' will be studied in detail for anisotropic traps. However, it will turn out that in traps there 
exists a second low-energy regime, which for isotropic traps is defined by E/ ^ — > with E — (ljqL)'^ /Afi <|; E, where 
the quasi-particles are single-particle like excitations confined to a narrow layer around the surface of the condensate. 
This low-energy Hartree-Fock regime will be discussed in detail in section IV, also for anisotropic traps, together with 
the usual high-energy {E ^ /i) Hartree-Fock regime. 



III. QUASI-PARTICLE DYNAMICS IN THE HYDRODYNAMIC REGIME 

A. Hydrodynamic Hamiltonian 

If there exist limiting trajectories for different initial conditions there should exist limiting dynamics described by 
some limiting Hamilton-function. Inside the condensate the Bogoliubov Hamiltonian can be written as 

H{p, x) = y/ek^„{p){ek^n{p) + ^2K{x)) , (27) 

where Ckmip) — p^/2m. For small energies K{x) is much bigger than ekin everywhere except in a small region near 
the boundary. This suggests that the approximant of the Hamilton function ( |27| ) can be obtained by neglecting the 
kinetic energy square 

HhydiP,x) = ^2ekrn{p)K{x) (28) 

for describing the motion in the hydrodynamical regime. This approximate Hamiltonian is in accordance with the 
bulk case, when K{x) = fi should be taken in ( p8| ) in order to obtain the linear phonon spectra from the Bogoliubov 
dispersion relation. 

This Hamiltonian is meaningful only inside the condensate and only near the boundary of the condensate the 
full Hamiltonian (^) differs from this approximate one. On the Thomas-Fermi surface the full Hamilton function 
gives definite values for the Cartesian momenta, whereas according to H^yd they become infinite. Following the 
trajectories of H^yd m the isotropic case the angular momentum conservation requires that the tangential component 
of the momentum remains finite even though the absolute value of the momentum diverges hke Therefore 
each trajectory hits the boundary orthogonally and is reflected back orthogonally without change in the tagential 
component of the momentum. As this local rule is independent of the global symmetry of the trap potential it must 
hold also in the anisotropic case. 

The Hamiltonian ( |2^ ) has some further unusual features. The first observation is that it is not of the usual form 
but is a homogeneous first order function of the momenta. The strong consequence is that with the same initial value 
x{t = 0) and with the same direction of the initial momenta the orbit x{t) is the same independently of the energy. 
Secondly, a constraint follows from the canonical equations of motion, namely 
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mx • X = fi — U{x) , 



(29) 



relating the velocities and the coordinates. Thus one cannot choose the initial point and the velocity independently. 
Furthermore, due to this constraint one cannot express the three velocities in terms of the momenta, i.e., one cannot 
do the inverse Legendre transformation in the usual way to derive the Lagrangian. From ( p9| ) it is clearly seen that 
despite of the divergence of momenta on the boundary of the condensate the velocities even tend to zero here. 



B. Isotropic case 

In the isotropic trap case (wo ) the Poisson-bracket of Hhyd{p,x) (See Eq.(p8|)) and the angular 

momentum vector L is zero, which means that any components of L is a conserved quantity. Let us choose our 
coordinate system in such a way that the z axis is parallel with L. In such a frame = Ly = 0, which shows that the 
motion in the phase space stays on the hypersurface z = pz = 0. By this choice of the coordinates one can eliminate 
one degrees of freedom from the Hamiltonian (Eq) , which has then the form 



HnydiP, {Pl + Pi) (1 - . (30) 



where ttf denotes the radial size of the condensate, the Thomas Fermi radius rxp =^ \ Let us now consider 

the transformation 

rTF 



ih + h) 

ih + h) 

ih + h) 



h cos 4>i + h cos 02) , 
Ii sin 01 - h sin 02), 

sin 01 + sin 02), 



?'Tf(1 - cos (01 + 02)) 



Py^ 77"^~^"~^7XT^rTr(*^°^'^i "*=°^'^2), (31) 

rTF(l-COS(0i+02)) 

with positive /i and l2- It is straightforward to check that the Poisson-brackets between /i,/2,0i,02 are canonical, 
thus, the transformation (^ij) is a canonical transformation. Inserting (|3|) into (^0|) one gets 



E = H{h,l2,^l,<l>2)=LOoV^hl2, (32) 



i.e., Ii and I2 are the action and 0i and 02 the angle coordinates of the Hamiltonian (30). Similarly to the harmonic 
oscillator case this Hamiltonian is a homogeneous first order function of the action coordinates. 
The Hamilton equations in the new coordinates are 

72 = 02 = i^i,/Il = ni. (33) 

V2 V -'2 

Using the above transformation, it is easy to show that the angular mumentum is 

Lz ^ xpy - ypx = h- I2. (34) 

A nice geometrical meaning for x{t) can be given. Let us consider a circle of radius 6, in which a smaller circle of 
radius a rolls. The motion of a point on the perimeter of the smaller circle in cartesian coordinates is described by 
the equations 

X = {b — a) cos 01 + a cos 02 , 

y ~ (b — a) sin 0i — a sin 02 , (35) 

where 0i and 02 are linear functions of the time, see Fig.|^. Due to the perfect rolling condition the angular velocities 
are not independent: 
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= {b - a)(^i - a(f>2. (36) 
Comparing the parametric form of the hypocycloid ( |35|) with ( ^l| ) it is obvious that x{t) fulfills ( ^5|) and the constraint 



(|36|), if = ttf, a = ?'tf/2/(/i + /2), and if 02 are chosen as in (33). 
The radial distance from the origin can be expressed by 



r = ^x^ + y2 ^ -p^Jlf + q + 2/1/2 cos (0i + ^2). (37) 

+ -12 

It is obvious that it is periodic in (0i + (^2), its period can be calculated from 2tt = (f2i + D,2)Tr, which yields 

Tr^- , ^ — (38) 
^0 ^2E^ + (L,^o)2 

The Hamiltonian ( ^0[ ) can be written in polar coordinates r, (/) as well. is a cyclic variable, its conjugate momentum 
^0 = \Lz\ is a conserved quantity. However, the momentum pr — {xpx + yPy)/r conjugated to r is not conserved. To 
express the Hamiltonian ( ^ ) in the action variables Ir , and /^ let us use the fact that 27r/r ^ § Pr dr and that during 
one period of the radial motion 0i + 02 changes by 27r. Using the above formulas one gets 

Ir min(/i,/2), 

/0 = 1/1-/2 1, (39) 

which leads by ( p^ ) to 

S = Hhydilr.U) - Uo^2{lr+I^)lr- (40) 

If one quantizes semiclassically the Hamiltonian ( ^ ) one should take into account that in the radial direction there 
are two turning points, thus, Ir should be replaced by ?i(n + 1/2) , and by the usual procedure for spherically symmetric 
problems /^ by h{l + 1/2) {I and n are non-negative integers). The semiclassical quantization leads by the above 
replacement rules to 



En,i = hujQ\j2n^ + 2nl + 3n + ^ + 1, (41) 

which is almost that of the result of Stringari [|l^, except the 1 under the square-root, and agrees with that of the 
more elaborate semiclassical quantization in the hydrodynamical limit . 



C. Anisotropic case with cylindrical symmetry 

The case of a trap with axial or cylindrical symmetry is the experimentally most relevant one. In Poincare cuts 
of the full dynamics we have seen regular behaviour for small energies. Therefore one can expect that the classical 
motion given by the approximate Hamiltonian is fully integrable. To show this let us introduce new coordinates, 
namely the cylindrical elliptical coordinates i^, r\ given by 

p = a^fWTW^) , z = a^Ti (42) 

which are orthogonal coordinates. Surfaces of constant ^ are confocal ellipsoids with foci at a distance a in p direction, 
surfaces of constant rj are confocal hyperboloids with the same foci. For a, the parameter of the transformation, we 
take the foci of the Thomas-Fermi ellipsoid, a = e{2p/mLo^Y/'^ for ujz > wo- For loq > ujz one has to change -f 1 
to — 1 and take a — e(2/i/TOw^)^/^. In the following we consider only the first case (p2[), in the second case the 
analysis proceeds similarly. ^ can take any value in the range [0, (1/e^ — 1)^^^]. The limiting case ^ = (1/e^ — 1)^/^ 
describes the Thomas- Fermi ellipsoid. 77 can be in the range [—1,1]. Making the point transformation from cylindrical 
to cylindrical elliptical coordinates the momenta transform as 

Pp = ^^2^v/(^' + - - Wv) , 

P.^-7I^ii^' + ^)m + i^~v')^Pn) ■ (43) 
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The Hamiltonian (E8^ in cylindrical elliptical coordinates is 



„2 + 

hyd ^2 + ^2 

X (ie + M + (1 - ^')pI + - ^)pi) ■ (44) 

Taking the energy E and = Lz as constants one can write down the Hamilton- Jacobi equation for ^ and 77, which 
is separable in these coordinates. Thus the problem is fully integrable. Introducing a separation constant B > the 
two separated Hamilton-Jacobi equations are 

e + 1 ujI 1-62(^2 + 1) 

Combining these two equations one gets for the separation constant B the phase-space function 

1 1 



B = 



g2 ^2 _|_ jj2 



+ (1 - e^{l mi v')pI + ij^ ^^)pI 



= -^[pl+Pl + il- ^^)pI] - (xPx + ypy + zpzf . (46) 

in elliptical and cartesian coordinates respectively. This is the third independent constant of motion in addition to the 
energy E and Lz ■ This can be checked directly, using the equations of motion for the time derivatives of B . Similarly 
to the isotropic case conservation of E and B require that trajectories hit the boundary orthogonally, because the 
momenta there diverge. In the isotropic limit a —f Q the elliptical coordinates become singular, and therefore it is 
more instructive to see this limit in cartesian coordinates. In this limit cr/e is the Thomas- Fermi radius, and B has 
the simple meaning 

B = ^ + L\ (47) 
^0 

The existence of the three independent constants of motion E, Lz and B explains the integrable motion generated by 
Hhyd and therefore the almost integrable situation found numerically in the motion generated by the total Hamiltonian 
( P7[ ) in the small energy and small angular momentum region. We notice that two kinds of trajectories can occur in 
this regime. From ( [45|) we can determine the turning points in ^ and 77. In ^-direction all the trajectories reach the 
Thomas-Fermi surface and are reflected back there. If the condition 

B> B* = + Ll , (48) 
^0 

is satisfied, there is an inner turning point in ^-direction and 77 takes a range [—rj^^^^, ?7max]- These trajectories correspond 
to the hypocycloids of the isotropic case, as an example see FigJJa. For B < B* however there are further turning 
points in 7;-direction, the motion being confined between two hyperbolas with 1^ values extending to zero, which can be 
seen in Fig.^. This kind of trajectory only occurs in the anisotropic system. B = B* is the separatrix between these 
two types of motion. As usual this separatrix is structurally unstable against small nonintegrable perturbations of 
the integrable motion in the hydrodynamic limit. It plays a crucial role for the appearance of chaos in the Bogoliubov 
Hamiltonian as the energy is increased from values very small compared to /i, because it is destroyed and replaced 
by a chaotic separatrix layer, which is very narrow at first, but grows in width as the energy is increased. In Fig.|l]c 
two regular islands corresponding to the two kinds of trajectories and the chaotic separatrix layer between them can 
be seen. 
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D. Completely anisotropic case 



The analysis of the preceding section can be generahzed to the case of a completely anisotropic harmonic trap. The 
formulas become rather lengthy and we just indicate the essential steps. 
The trap potential is written in the form 



Uix) 



X 



with 



52 



(and cyclic) , 



(49) 



(50) 



We may assume a > b > c without restriction of generality. Then new elliptic coordinates ^, 77, C are introduced via 

(and cyclic) 



X — ±A 



'(a2+0(a2 + 7y)(a2+C) 



after which the potential reads 



(a2-62)(„2_c2) 



c/(C, 77,0-^1 + ^1)- 



(51) 



(52) 



The range of ^, 77, C is > ^ ^^2 > rj > —iP' > ^ > —o?. The old canonical momenta Px, Py, Pz are given in terms 
of the new ones by 



Px 



'(a2+g)(a2+r7)(a2 + C) 
(a2-62)(a2_c2) 



2p^ 

+2pc 



(&2 + 0(c' + 6 



(e- 


-m- 





(62^ 


f 77)(c2 


+ v) 


(77 


-0(^7- 


-0 




C)(cM 


-0" 


(C- 


-m- 





(53) 



(and cychc) . 

Then the Hamiltonian in the hydrodynamic limit can be written in terms of the new canonical coordinates and 
momenta. The Hamilton- Jacobi equation can be written in the nearly separated form 



= (77 - C) 



(cyclic) 



(54) 



4^ C 

In fact complete separation is achieved, because this equation is satisfied only by putting the angular brackets equal 
to A + BS, (and cyclic), 



(«'+0(6'+0(c'+0Pf + 



4^l e 



(55) 



(and cyclic) 

where A and B are two separation constants, which are the same for all three equations related by cyclic permutation. 
From these three equations A and B can be eliminated by multiplying the first with (77 — C), the second with (C ~ 0, 
the third with (^ — 77) and adding them. This gives, of course, back eq. (p^, which defines E = ifhyd in terms of 
the canonical variables. However, solving the three equations instead for A by eliminating B and E, and then for B 
eliminating A and E we obtain two new conserved phase-space functions. Translated back to Cartesian coordinates 
these read 



{[(6^ + c2)(a;2 - a^) + a^y^ + z^)]pl + 2a^yzpyp^] + (cyclic) 



and 



B = -{{x^ - a^)pl + 2yzpyp^] + (cychc) 



(56) 
(57) 

The conserved function i? is a simple generalization of the conservation law we already found in the case with 
cylindrical symmetry, whereas A corresponds to L\ . By a straightforward but lengthy calculation it can be checked 
that the Poisson brackets {iJhyd,^}, {ffhyd,-B}, {^j -B} all vanish. Therefore, the dynamics governed by i/hyd is still 
completely integrable even in this completely anisotropic case. 
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IV. THE HARTREE-FOCK DYNAMICS 



Another limiting case of the Bogohubov description of quasiparticles (H) consists in neglecting the hole-component 
Vijx) in the field operator ifix). The remaining component Uj{x) is then described by the Hartree-Fock Hamiltonian 
(|lO|). The interaction between particles is taken into account by the potential K{x), describing the mean interaction of 
one particle with all the other particles. Restricting ourselves to T = all those other particles are in the condensate. 
In the homogenous systems this approach simply results in a shift of the dispersion relation of noninteracting particles 
by the chemical potential ^. For spatially homogeneous Bose condensates and also Bose condensates in traps such a 
description can be applied for energies larger than the mean interaction energy given by /i. However, in traps there 
is even a regime for energies smaller than /i where the Hartree-Fock approximation applies [ pl| , namely in the case 
when the kinetic energy Ckinip) is large compared to the local mean interaction energy K{x). This can be satisfied in 
a layer around the surface of the Bose condensate where K{x) is very small. 

Using the Thomas-Fermi approximation for the wave function (H) the Bogoliubov Hamiltonian and the Hartree- 
Fock Hamiltonian coincide outside the condensate. Inside the condensate, if the kinetic energy ekin{p) is much larger 
than the potential term K{x), an expansion of the Bogoliubov Hamiltonian ( p7[ ) to first order in K{x) just gives the 
Hartree-Fock Hamiltonian 

H^, = ^ + \U{x)-^i\, (58) 

which is therefore valid, for tkin >> K{x), inside and outside the condensate. 

We now want to investigate the classical dynamics of this Hartree-Fock Hamiltonian. The isotropic problem is 
completely integrable again. As constants of motion we can take the energy, the modulus and the z-component of 
the angular momentum. We immediately turn to the classical dynamics of the anisotropic, but axially symmetric 
case in the trap potential ( p5| ) and consider it as a system with two degrees of freedom. The conserved angular 
momentum around the symmetry axes Lz enters only as a parameter. Again we investigate the dynamics by Poincare 
cuts, now taken on the Thomas-Fermi surface ^ = — 1)^/^ and parametrized by the second elliptical coordinate 
T] and its conjugate momentum p^. For energies much larger than the chemical potential the interaction with the 
condensate is only a small perturbation to the integrable motion in the harmonic trap and we observe quasi-regular 
behaviour. In this limit the Bogoliubov description of quasi-particles reduces to the Hartree-Fock description, the 
condition ekin{p) ^ K{x) being fulfilled for all trajectories, and the classical motions generated by both Hamiltonians 
are essentially the same. Trajectories not entering the condensate are even identical, since here the two descriptions 
fully coincide. 

For energies in the approximate range 10 > E / ^ > 0.1 we observe a mixed phase space again (see Fig.^). A regular 
island around the periodic orbit z = = pz is surrounded by a chaotic sea. For E > fi the structure in phase space 
is similar as for the Bogoliubov dynamics, but differing in detail. From this we can conclude that the stochasticity 
observed for the Bogoliubov Hamiltonian in Fig.|l| is not a consequence of the coupled two component structure of the 
underlying semiclassical description, but is simply caused by the anisotropy of the external potential. 

For energies much smaller than the chemical potential E < 0.1/i we find regular behaviour again, see Fig.^D. 
Particles are confined to the sharp potential channel near the Thomas-Fermi surface. The width of this potential 
channel scales as E/^. Roughly the particles spend the same time outside and inside the condensate. We look at 
the problem in elliptical coordinates (^), and choose for concreteness again the case lOz > ujq. The oscillations in ^ 
orthogonal to the Thomas-Fermi surface ^ = (1/e^ — 1)^/^ are much faster than the oscillations in 77-direction along 
the channel. This suggests to make an adiabatic approximation in which the action- integral = {2tt)^^ § P(,d£, over 
one full cycle in ^ at fixed 77, emerges as an adiabatic constant for the motion. Evaluating this adiabatic invariant 
for £'//i <C 1 we get as a function of 77, p,, 

This new adiabatically conserved quantity which emerges in the low-energy limit of the Hartree-Fock dynamics, is 
the cause of integrability in that limit. 

Solving this equation for the energy we get the Hamiltonian of the slow 7/-dynamics, valid for low energies. From 
( p9| ) we see that the turning points in rj are independent of the energy if we keep Lz/E'^ and I^/E^^'^ constant. In 
Fig.^ we compare a Poincare section in rj and of the dynamics of the Hartree-Fock Hamiltonian with trajectories 
of the slow //-dynamics for different values of /j. Both curves agree very well. For smaller energies E < 0.01/i no 
difference between both curves can be noticed. 

Now we have to ask ourselves, which of the trajectories displayed in Fig.^ are indeed good approximations to 
trajectories described by the Bogoliubov Hamiltonian. Let us look first at the isotropic case, where the motion 
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separates in radial and angular motion. The kinetic energy in the angular degree of freedom roughly is L^/2mr^ « 
{uJoL)^ /Ajj,. Since for low energies r deviates only very little from the Thomas-Fermi radius this rotational energy 
is almost conserved. The remaining energy is stored in the radial degree of freedom and only this energy can be 
transformed to potential energy. So the condition that the Hartree-Fock dynamics and the Bogoliubov dynamics 
agree is in this case 

(60) 

For the anistropic case we can formulate an analogous criterion: Only if most of the energy is kinetic energy of 
the motion parallel to the boundary, which cannot be transformed into potential energy, the approximation of the 
Bogoliubov dynamics by the Hartree-Fock dynamics works well. For = no motion takes place orthogonal to the 
surface, being constant, and all the energy is stored in motion parallel to the surface. This corresponds to the outer 
orbit forming the boundary of the cut in Fig.^. The maximal value of for fixed Lz is given by setting 77 = p,, = 
in (^9|). This trajectory corresponds to the origin of Figj^, where motion takes place only in ^ and 0-directions. 
Generally the two dynamics differ in this case, unless most of the energy is stored in angular motion of the cyclic 
variable (f) which is also motion along the Thomas-Fermi surface. The maximal value of /{ for a given energy is found 
by neglecting both angular motions in 77 and in (p in (|59|), by putting p^, equal to zero there. As a condition that 
only a small fraction of energy is stored in the motion orthogonal to the surface and hence that both motions from 
the Bogoliubov and from the Hartree-Fock Hamiltonian agree, we can thus state: 

/5«/:'-- = -^f-)'^'. (61) 



SvrijJo V ji 

In Fig.O we compare Poincare cuts of the Bogoliubov dynamics with the one-degree of freedom motion obtained 
from (p9h, representing the integrable Hartree-Fock dynamics for small energies. We see that indeed both dynamics 
agree well for small values of near the boundary of the cut. For chosen large, see Figj^a, even for values of 
close to the maximal one, both dynamics in 77 agree qualitatively. However the different behaviour in the variable ^ 
orthogonal to the surface can, of course, not be seen in this cut at constant ^. For smaller values of L^, see Fig.^, 
we can distinguish two regions. Near the boundary, for small we see the Hartree-Fock limit of the Bogoliubov 
dynamics, where both dynamics agree. The inner region corresponds to the hydrodynamic limit of the Bogoliubov 
dynamics and cannot be compared with the Hartree-Fock dynamics. The two kinds of closed tori visible here are the 
two kinds of hydrodynamic trajectories discussed at the end of section HI.C. 



V. CONCLUSIONS 



The quasi-particle excitations are the basic constituents of the dynamical and thermodynamical properties of Bose 
condensates. In the present paper we have investigated their dynamics for Bose condensates of atomic gases in traps 
in the classical limit. The two limiting types of excitations, collective modes and quasi-particle excitations consisting 
essentially of single atoms moving in a mean field correspond, in the classical limit, to particles and anti-particles of 
zero mass, moving 'relativistically' with the speed of sound, and to single atoms moving in the potential created by 
the trap and the Hartree-Fock potential energy of all other atoms. In spatially homogeneous (untrapped) condensates 
these two types of excitation strongly differ in energy E, the collective modes occurring at i? ^ /x, the single-particle 
modes at E ^ fj,. In the trapped condensates both types of excitations coexist, at least classically, at small energies 
E <^ fj,, and are instead spatially separated. The collective modes live inside the condensate, the single-particle modes 
at small energies in a narrow layer at the border. 

One principal result we have obtained here is that the classical dynamics of both, the collective modes and the 
single-particle modes, become integrable in the limit E/fi ^ 1. This has important consequences for the quantum 
dynamics as well: the integrability can be used there to separate the Schrodinger equation and to obtain not only the 
low- lying levels of the collective modes ||l^,|l^, but also of the single-particle modes. After quantization an energy 
gap reappears separating the collective modes with typical energies hujo and the single-particle modes whose lowest 
levels have energies of the order (Tiwo)^^'^/^^^^ due to their close confinement in normal direction to the surface of the 
condensate. However, this energy difference is much smaller than, and has a different origin as the energy difference 
between both types of modes in homogeneous systems. 

Another principal result obtained here is the nonintegrabilty of the classical dynamics of the quasi-particle excitations 
at intermediate energies E :^ fx. This applies to both, the full Bogoliubov dynamics and the limiting Hartree-Fock 
dynamics approximating it wherever the kinetic energy is large compared to the local mean interaction energy. Again 
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this nonintegr ability has a direct consequence also for the quantum dynamics, because it implies avoided crossings 
between quasi-particle levels as functions of the dimensionless interaction strength Nga/do with do = y^h/rriLLiQ, if the 
energy and fj, are comparable. Such avoided crossings have indeed been seen in numerically generated plots . 

Our results not only explain these avoided crossings, they also open the door to an intreaguing wider perspective, 
quantum chaos of the quasi-particle dynamics in the Bose condensates of atoms in anisotropic traps. 
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FIG. 1. Poincare sections of the dynamics of the Bogoliubov Hamiltonian (^]) in cylindrical coordinates for the different 
energies (from top to bottom) E/jj, = 40 (a), 1 (b) and 0.02 (c). The cut is taken at z — and diplayed in the variables p,Pp 
in units of (2/i/ma;o)^'''^, (2m/i)^''^, respectively. The anisotropy is chosen as bj^jijjQ = VS, the angular momentum was fixed 
as ujqL^/E — 0.2. 

FIG. 2. Trajectories in coordinate space of the Bogoliubov dynamics of (Gil) starting from the same point and in the same 
direction for different energies E/fi — 0.1 (dashed line), 0.01 (solid) and 10~^dotted). 

FIG. 3. Coordinates for the hypocycloid (paj) 
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FIG. 4. Trajectory of the hydrodynamic Hamiltonian ( |28| ) for B > -B*(a) and for B < B*{h). z,p are plotted in units of 



FIG. 5. Poincare section of the dynamics of the Hartree-Fock Hamiltonian (jSg) in elliptical coordinates ( |42[ ) for the energy 
E/^ — 1(a) and E/^ = .06(b). The cut is taken on the Thomas-Fermi surface in the variables 77,^,7 for uJz/luq = \/8. The 
angular momentum is given by (ujoL^)^ /2fj.E = 1. Solid lines in (b) are trajectories of the Hamiltonian in {rj,Pr,) following from 
(§)• 

FIG. 6. Poincare sections of the dynamics of the Bogoliubov Hamiltonian (^l|) in elliptical coordinates (^) for the energies 
E/fi = 0.1 with (ujoLz)^ /2fj,E = 1(a) and E/fi — 0.01 with (uJoLz)^ /2fiE = .02(b) on the Thomas-Fermi surface in r],Pn for 
i^z/i^o = VS- Solid lines are trajectories of the Hamiltonian in {r],p,^) following from (^j]). 
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